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Abstract. We review mechanisms of low-temperature electronic transport 
through a quantum dot weakly coupled to two conducting leads. Transport in this 
case is dominated by electron-electron interaction. At temperatures moderately 
lower than the charging energy of the dot, the linear conductance is suppressed 
by the Coulomb blockade. Upon further lowering of the temperature, however, 
the conductance may start to increase again due to the Kondo effect. We 
concentrate on lateral quantum dot systems and discuss the conductance in a 
broad temperature range, which includes the Kondo regime. 



PACS numbers: 72.15.Qm 73.23.-b 73.23.Hk 73.63.Kv 

1. Introduction 

In quantum dot devices |T] a small droplet of electron liquid is confined in a finite region 
of space. The droplet can be attached by tunneling junctions to massive electrodes 
to allow electronic transport across the system. The conductance of such a device 
is determined by the number of electrons on the dot N, which in turn is controlled 
by varying the potential on the gate - an auxiliary electrode capacitively coupled to 
the dot pp. At sufficiently low temperatures the number of electrons N is an integer 
at almost any gate voltage Vg. Exceptions are narrow intervals of Vg in which an 
addition of a single electron to the dot does not change significantly the electrostatic 
energy of the system. Such a degeneracy between different charge states of the dot 
allows for an activationless electron transfer through it, whereas for all other values of 
Vg the activation energy for the conductance G across the dot is finite. The resulting 
oscillatory dependence G{Vg) is the hallmark of the Coulomb blockade phenomenon 
The contrast between the low- and high-conductance regions (Coulomb blockade 
valleys and peaks, respectively) gets sharper at lower temperatures. This pattern 
of G{Vg,T) dependence is observed down to the lowest attainable temperatures in 
experiments on tunneling through small metallic islands [2]- However, small quantum 
dots formed in GaAs heterostructures exhibit drastically different behavior |3| : in some 
Coulomb blockade valleys the dependence G{T) is not monotonic and has a minimum 
at a finite temperature. This minimum is similar in origin ^ to the well-known non- 
monotonic temperature dependence of the resistivity of a metal containing magnetic 
impurities - the Kondo effect. 

In this paper we review the theory of the Kondo effect in quantum dots, 
concentrating on the so-called lateral quantum dot systems formed by 

gate depletion of a two-dimensional electron gas at the interface between two 
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semiconductors. These devices offer the highest degree of tunabiUty, yet allow for 
relatively simple theoretical treatment. At the same time, many of the results 
presented below are directly applicable to other systems as well, including vertical 
quantum dots El El j Coulomb-blockaded carbon nanotubes [3 (HI , single- molecule 
transistors |10| . and stand-alone magnetic atoms on metallic surfaces 

2. Model of a lateral quantum dot system 

The Hamiltonian of interacting electrons confined to a quantum dot has the following 
general form, 

^dot = E E ^^Asd.s + ^ E E K,Md\sd],d,,d,^. (2.1) 

s ij ss' ijkl 

Here an operator d\^ creates an electron with spin s in the orbital state (j)i{r); h^ - — h*- 
is an Hermitian matrix describing the single-particle part of the Hamiltonian. The 
matrix elements hijki depend on the potential U (r—r') of electron-electron interaction, 

h,,ki = I drdr'c^*{r)c^*{r')U{r - r')Mr>i (r). (2.2) 

The Hamiltonian (|2.1|l can be simplified further provided that the quasiparticle 
spectrum is not degenerate near the Fermi level, that the Fermi-liquid theory is 
applicable to the description of the dot, and that the dot is in the metallic conduction 
regime. The first of these conditions is satisfied if the dot has no spatial symmetries, 
which implies also that motion of quasiparticles within the dot is chaotic. 

The second condition is met if the electron-electron interaction within the dot is 
not too strong, i.e. the gas parameter Ts is small, 

r, = {kpao)-'^ < 1, ao = nh^/e^m* (2.3) 

Here kp is the Fermi wave vector, is the effective Bohr radius, k, is the dielectric 
constant of the material, and m* is the quasiparticle effective mass. 

The third condition requires the ratio of the Thouless energy Et to the mean 
single-particle level spacing 5E to be large jl2| . 

g = Et/5E > 1. (2.4) 

For a ballistic two-dimensional dot of linear size L the Thouless energy Et is of the 
order of Hvp/L, whereas the level spacing can be estimated as 

SE - hvphp/N - tf/m*L^. (2.5) 

Here is the Fermi velocity and N ~ {kpL)^ is the number of electrons in the dot. 
Therefore, 

g - kpL , 

so that having a large number of electrons TV ^ 1 in the dot guarantees that the 
condition (|2.4|) is satisfied. 

Under the conditions (|2.3|) , (|2.4|) the Random Matrix Theory (see QSl Ql for 
a review), is a good starting point for description of non- interacting quasiparticles 
within the energy strip of width Et about the Fermi level ^2 . The matrix elements 
hij in Eq. H2.1|l belong to a Gaussian ensemble Since the matrix elements do 

not depend on spin, each eigenvalue e„ of the matrix hij represents a spin-degenerate 
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energy level. The spacings e„+i — e„ between consecutive levels obey the Wigner-Dyson 
statistics Jl]; the mean level spacing e„+i — e„ — 5E. 

We discuss now the second term in the Hamiltonian l|2.1(l . which describes 
electron-electron interaction. It turns out |15l 1161 [T7] that the vast majority of the 
matrix elements hijki are small. Indeed, in the lowest order in \/g <^ 1, the wave 
functions 4>i{r) are Gaussian random variables statistically independent of each other 
and of the corresponding energy levels |18) : 



</,*(r)^^.(r') = j^S,,F{\r - r'|), F(r) ~ (exp(ifc • r))FS 



(2.6) 



Here (. . .)fs stands for the averaging over the Fermi surface |fc| = kp. In two 
dimensions, the function F{r) decreases with r as F cx (kpr)^^^^ at kpr 3> 1, and 
saturates to ^ 1 at kpr <C 1. After averaging with the help of Eq. (|2.6() . the matrix 
elements H2.2II take the formf 



hijki = "^EcSiiSjk + EsSikSji. 

We substitute this expression into the Hamiltonian (|2.1|l , and rearrange the sum over 
the spin indices with the help of the identity 

2 5„ , (5,',' — s' + <T„ „' • (T„ „' , (2.7) 

^1^2 12 11 ^2 2 11 ^2 2' ^ ' 

where <t = (tr^, cr^, cr^) are the Pauli matrices. This results in a remarkably simple 
form [Hini 

iJint = EcN^ - EsS^ (2.8) 
of the interaction part of the Hamiltonian of the dot. Here 

A^=>'dLd_, S = ^dL^d„,, (2.9) 



ns 



,1 



are the operators of the total number of electrons in the dot and of the dot's spin, 
respectively. 
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Figure 1. Equivalent circuit for a quantum dot connected to two leads by 
tunneling junctions and coupled via a capacitor (with capacitance Cg) to the 
gate. The total capacitance of the dot C = Cl + Cr + Cg. 

The first term in Eq. (|2.8|) represents the electrostatic energy. In the conventional 
equivalent circuit picture, see Fig. ^ the charging energy Ec is related to the total 
capacitance C of the dot, Ec = e^/2C. For a mesoscopic [kpL ^ 1) conductor, the 

t For simplicity we assumed here that (f>i{r)if>j {r') = 0, which corresponds to broken time-reversal 
symmetry. See 1171 for discussion of the general case. 
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charging energy is large compared to the mean level spacing 5E. Indeed, using the 
estimates C ^ kL and (|2.5|) . we find 



Ec/6E ^ L/ao rWN. (2.10) 

Except an exotic case of an extremely weak interaction, this ratio is large for iV ^ 1; 
for the smallest quantum dots formed in GaAs heterostructures, Ec/SE ~ 10 

The second term in Eq. (|2.8|l describes the intra-dot exchange interaction, with 
the exchange energy Eg given by 

Es = J drdr'U{r-r')F^{\r-r'\) (2.11) 

In the case of a long-range interaction the potential U here should properly account for 
the screening ^7]. For ^ 1 the exchange energy can be estimated with logarithmic 
accuracy by substituting U{r) — [e^ / nr)9{ao — r) into Eq. H2.11|l (here we took into 
account that the screening length in two dimensions coincides with the Bohr radius 
oq), which yields 

Es ^rs\iv{l/rs)5E -^SE. (2.12) 

The estimate (|2.12|l is valid only for ^ 1. However, the ratio Es/SE remains small 
for experimentally relevant| value ~ 1 as long as the Stoner criterion for the absence 
of itinerant magnetism is satisfied. This guarantees the absence of a macroscopic 
(proportional to N) magnetization of a dot in the ground state [16j . 

Obviously, the interaction part of the Hamiltonian, Eq. I|2.8|l . is invariant with 
respect to a change of the basis of single-particle states (jiiir). Picking up the basis in 
which the first term in H2.1|l is diagonal, we arrive at the universal Hamiltonian |16[ll7j . 

iJdot =Y.^rAsdns+Ec{N~Nof - EsS\ (2.13) 

ns 

We included in Eq. (|2.13|l the effect of the capacitive coupling to the gate electrode: the 
dimensionless parameter is proportional to the gate voltage, = CgVg/e, where 
Cg is the capacitance between the dot and the gate, see Fig.^ The relative magnitude 
of various terms not included in (|2.13() , as well as that of mesoscopic fluctuations of 
the coupling constants Ec and Es, is of the order oil/g ^ N^'^^^ ^ 1. 

As discussed above, in this limit the energy scales involved in 12.13|l form a well- 
defined hierarchy 

Es<t:SE<t: Ec. (2.14) 

If all the single-particle energy levels e„ were equidistant, then the spin S of an even-A^ 
state would be zero, while an odd-A state would have S — 1/2. However, the level 
spacings are random. If the spacing between the highest occupied level and the lowest 
unoccupied one is accidentally small, than the gain in the exchange energy, associated 
with the formation of a higher-spin state, may be sufficient to overcome the loss of 
the kinetic energy (cf. the Hund's rule in quantum mechanics). For Es «C 6E such 
deviations from the simple even-odd periodicity are rare [161 120L |2] ■ This is why the 
last term in H2.13|) is often neglected. Equation (|2.13|) then reduces to the Hamiltonian 
of the Constant Interaction Model, widely used in the analysis of experimental data ^ . 

f For GaAs (m* 0.07me , k 13) the effective Bohr radius ao ^ 10 nm, whereas a typical density 
of the two-dimensional electron gas, n ~ 10^^ cm~^ :3 , corresponds to kp - 
This gives kpao ~ 1. 
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Electron transport through the dot occurs via two dot-lead junctions. In a typical 
geometry, the potential forming a lateral quantum dot varies smoothly on the scale of 
the Fermi wavelength, see Fig. |21 Hence, the point contacts connecting the quantum 
dot to the leads act essentially as electronic waveguides. Potentials on the gates control 
the waveguide width, and, therefore, the number of electronic modes the waveguide 
support: by making the waveguide narrower one pinches the propagating modes off 
one-by-one. Each such mode contributes 2e^//i to the conductance of a contact. 
The Coulomb blockade develops when the conductances of the contacts are small 
compared to 2e^//i, i.e. when the very last propagating mode approaches its pinch- 
off 22, 23|. Accordingly, in the Coulomb blockade regime each dot-lead junction in a 
lateral quantum dot system supports only a single electronic mode |24) . 




Figure 2. The confining potential forming a lateral quantum dot varies smoothly 
on the scale of the de Broglie wavelength at the Fermi energy. Hence, the dot-lead 
junctions act essentially as electronic waveguides with a well-defined number of 
propagating modes. 

As discussed below, for Ec ^ 5E the characteristic energy scale relevant to the 
Kondo effect, the Kondo temperature Tk, is small compared to the mean level spacing: 
Tk <C 5E. This separation of the energy scales allows us to simplify the problem even 
further by assuming that the conductances of the dot-lead junctions are small. This 
assumption will not affect the properties of the system in the Kondo regime. At the 
same time, it justifies the use of the tunneling Hamiltonian for description of the 
coupling between the dot and the leads. The microscopic Hamiltonian of the system 
can then be written as a sum of three distinct terms, 

H = iJlcads + Hdot + ^^tunncling, (2-15) 

which describe free electrons in the leads, an isolated quantum dot, and tunneling 
between the dot and the leads, respectively. The second term in H2.15|l . the 
Hamiltonian of the dot -ffdot: is given by Eq. (|2.13|) . We treat the leads as reservoirs of 
free electrons with continuous spectra ^fc, characterized by constant density of states 
V, same for both leads. Moreover, since the typical energies lo < Ec of electrons 
participating in transport through a quantum dot in the Coulomb blockade regime 
are small compared to the Fermi energy of the electron gas in the leads, the spectra ^fc 
can be linearized near the Fermi level, = vpk; here k is measured from kp- With 
only one electronic mode per junction taken into account, the first and the third terms 
in Eq. H2.15|l have the form 

iJloads ^^^k^lks^aks' = -£.-k, (2.16) 

aks 

i^tunneling — ^an^aks^ns + H.C. (2-17) 

akns 
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Here tan are tunneling matrix elements (tunneling amplitudes) "connecting" the state 
n in the dot with the state k in the lead a {a = R,L for the right/left lead). The 
randomness of states n translates into the randomness of the tunneling amplitudes. 
Indeed, the amplitudes depend on the values of the electron wave functions at the 
points Ta of the contacts, t 

an Foi kpL ^ 1 tliG wavc functions are Gaussian 

random variables. Equation (|2.6(1 then results in 



iS„ia'„' = \tln\Sac.'6nn'- (2-18) 

Average values of the tunneling probabilities can be expressed via the conductances 
of the dot-lead junctions Ga, 



Here is the average rate for an electron to escape from a discrete level n in the dot 
into lead a. 

3. Rate equations and conductance across the dot 

At high temperatures, T 3> Ec^ charging energy is negligible compared to the thermal 
energy of electrons. Therefore the conductance of the device in this regime Goo is not 
affected by charging and, independently of the gate voltage, 

111 

7^ = 7^ + 7^- (3-1) 

Lfoo LfL LriJ 

Dependence on iVo develops at lower temperatures, 

SE<^T-^Ec. (3.2) 

The conductance is not suppressed only within narrow regions - Coulomb blockade 
peaks, i.e. when the gate voltage is tuned sufficiently close to one of the points of 
charge degeneracy, 

\N,~N*\<T/Ec; (3.3) 

here A^g is a half-integer number. 

We will demonstrate this now using the method of rate equations • In addition 
to constraint 13.2|l . we will assume that the inelastic electron relaxation rate within 
the dot is large compared to the escape rates Ta- In other words, transitions between 
discrete levels in the dot occur before the electron escapes to the leads§. Under this 
assumption the tunnelings across the two junctions can be treated independently of 
each other. Condition H3.3|l . on the other hand, allows us to take into account only two 
charge states of the dot which are almost degenerate in the vicinity of the Coulomb 
blockade peak. For A'q close to Nq these are the states with A^ = A'q ±1/2 electrons on 
the dot. Hereafter we denote these states as |A^) and |A^-|- 1). Finally, condition H3.2|l 
enables us to replace the discrete single-particle levels within the dot by a continuum 
with the density of states 1/SE. 

Applying the Fermi Golden Rule and using the described simplifications, we may 
write the current from the lead a into the dot as 

= ^ ^ \tan\'^5{S,k + cVa + En ~ e„ - Ejy+i) 
kns 

X {PNf{^k)[l - /(£«)] - PN+lf{en)[l - /(6)]} . 

§ Note that a finite inelastic relaxation rate requires inclusion of mechanisms beyond the model 12.131 . 
e.g., electron-phonon collisions. 
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Here /(i^) — [exp(w/T) + 1]~^ is the Fermi function, is the potential on the 
lead a, see Fig. ^ Ejq and i?7v+i are the electrostatic energies of the charge states 
|A^) and |7V + 1), and P/y and Pn+i are the probabilities to find the dot in these 
states. Replacing the summations over n and k by integrations over the corresponding 
continua, we find 

/a = — [PnF{Ei -Eo- eV^) - PN+iFiEo - E^ + eV^)] (3.4) 
e 

with F{uj) — w/[exp(a;/r) — f ]. In equilibrium, the current (|3.4() is zero by the detailed 
balance. When a finite current flows through the junction the probabilities Pjv and 
Pjv+i deviate from their equilibrium values. In the steady state, the currents across 
the two junctions satisfy 

I = Il = ~Ir- (3.5) 

Equations H3.4|l and (|3.5|l . supplemented by the obvious normalization condition 
Pjv + Pn+1 — 1, allow one to find P/v, Pw+i, and the current across the dot / in 
response to the applied bias V = Vl — Vr. The linear conductance across the dot 
G = limy^o dl/dV is then given by |^ 

Ec{N,-N*,)/T 

''-^°",nM2Ec{N,-N*)/T]- ^''^^ 

Here A'q — — (half- integer A^o) corresponds to the Coulomb blockade peak. In 
the Coulomb blockade valleys (Nq ^ Nq), conductance falls off exponentially with the 
decrease of temperature, and all the valleys behave exactly the same way. 



4. Activationless transport through a blockaded quantum dot 

According to the rate equations theory |25|. at low temperatures, T <^ Ec, conduction 
through the dot is exponentially suppressed in the Coulomb blockade valleys. This 
suppression occurs because the process of electron transport through the dot involves 
a real transition to the state in which the charge of the dot differs by e from 
the thermodynamically most probable value. The probability of such fluctuation is 
proportional to exp {—Ec\Nq — N^\/T), which explains the conductance suppression, 
see Eq. (|3.5f) . Going beyond the lowest-order perturbation theory in conductances 
Ga allows one to consider processes in which states of the dot with a "wrong" charge 
participate in the tunneling process as virtual states. The existence of these higher- 
order contributions to the tunneling conductance was envisioned first by Giaever and 
Zeller The first quantitative theory of this effect, however, was developed much 
later 

The leading contributions to the activationless transport, according to |2Z], are 
provided by the processes of inelastic and elastic co-tunneling. Unlike the sequential 
tunneling, in the co-tunneling mechanism, the events of electron tunneling from one of 
the leads into the dot, and tunneling from the dot to the other lead occur as a single 
quantum process. 



^.1. Inelastic co-tunneling 

In the inelastic co-tunneling mechanism, an electron tunnels from a lead into one of 
the vacant single-particle levels in the dot, while it is an electron occupying some 
other level that tunnels out of the dot, see Fig. |SIa). As a result, transfer of charge e 
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between the leads is accompanied by a simultaneous creation of an electron-hole pair 
in the dot. 




(b) 





Figure 3. Examples of the co-tunneling processes. 

(a) inelastic co-tunneling: transferring of an electron between the leads leaves 
behind an electron-hole pair in the dot; (b) elastic co-tunneling; (c) elastic co- 
tunncling with a flip of spin. 



Here we will estimate the contribution of the inelastic co-tunneling to the 
conductance deep in the Coulomb blockade valley, i.e. at almost integer Nq. Consider 
an electron that tunnels into the dot from lead L. If energy lo of the electron relative 
to the Fermi level is small compared to the charging energy, to Ec, then the energy 
of the virtual state involved in the co-tunneling process is close to Ec- The amplitude 
Ain of the inelastic transition via this virtual state to lead R is then given by 



A ■ — 



Rn' 



Ec 



(4.1) 



The initial state of this transition has an extra electron in the single-particle state k 
in lead L, while the final state has an extra electron in the state k' in lead R and an 
electron- hole pair in the dot (state n is occupied, state n' is empty). 

Given the energy of the initial state uj, the number of available final states can 
be estimated from the phase space argument, familiar from the calculation of the 
quasiparticle lifetime in the Fermi liquid theory 28 . For uj 3> SE this number is of 
the order of (uj/SE)'^. Since the typical value of lu is T, the inelastic co-tunneling 
contribution to the conductance can be estimated as 



Gi' 



T y 

SE 



(4.2) 



Using now Equations 12.1811 and (|2.19|) , we find |27] 

Gin ~ ~^ GlGji —— 

& \Ec , 

A comparison of Eq. 14.2|l with the result of the rate equations theory 1)3. 6|l shows that 
the inelastic co-tunneling takes over the thermally-activated hopping at moderately 
low temperatures 



T < Tu 



Er 



In 



(4.3) 



\ Gl + Gr , 

The smallest energy of the electron-hole pair is of the order of 5E. At 
temperatures below this threshold the inelastic co-tunneling contribution to the 
conductance becomes exponentially small. It turns out, however, that even at 
much higher temperatures this mechanism becomes less effective than the elastic co- 
tunneling. 
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4-2. Elastic co-tunneling 

In the process of elastic co-tunneling, transfer of charge between the leads is not 
accompanied by the creation of an electron-hole pair in the dot. In other words, 
occupation numbers of single-particle energy levels in the dot in the initial and final 
states of the co-tunneling process are exactly the same, see Fig. |2Ib). Close to the 
middle of the Coulomb blockade valley (at almost integer Nq) the average number of 
electrons on the dot, N k, Nq, is also an integer. Both an addition and a removal of 
a single electron cost Ec in electrostatic energy, see Eq. (|2.13|l . The amplitude of the 
elastic co-tunneling process in which an electron is transfered from lead L to lead R 
can then be written as 

A^,=Y^A^, A„=tl„t^„%l-fc^. (4.4) 

n ' ' 

The two contributions to the partial amplitude An are associated with virtual creation 
of either an electron if the level n is empty (e„ > 0), or of a hole if the level is occupied 
(e„ < 0); the relative sign difference between the two contributions originates in the 
fermionic commutation relations. 

With the help of Equations (|2.18|) and H2.19|) the average value of the elastic 
co-tunneling contribution to the conductance can be written as 

=2 



Since for Ec S> 5E the number of terms making significant contribution to the sum 
over n here is large, and since the sum is converging, one can replace the summation 
by an integral which results in P2| 

G7i^ -^GlGb.^. (4.5) 
Ec 

Comparison with Eq. 14. 2() shows that the elastic co-tunneling dominates the electron 
transport already at temperatures 



T<Tei^ y/EcSE , (4.6) 

which may exceed significantly the level spacing. Note, however, that mesoscopic 
fluctuations of Gei are strong ^29^, of the order of its average value (|4.5() . Thus, 
although Gez is always positive, see Eq. H4.6|l . the sample-specific value of Gei for a 
given gate voltage may vanish [M)\ . 



5. Effective low-energy Hamiltonian 

In the above discussion of the elastic co-tunneling we made a tacit assumption that all 
single-particle levels in the dot are either empty or doubly occupied. This, however, 
is not the case when the dot has a non-zero spin in the ground state. A dot with odd 
number of electrons, for example, would necessarily have a half-integer spin S. In the 
most important case of 5 = 1/2 the top- most occupied single-particle level is filled 
by a single electron and is spin-degenerate. This opens a possibility of a co-tunneling 
process in which a transfer of an electron between the leads is accompanied by a flip 
of electron's spin with simultaneous flip of the spin on the dot, see Fig. EJc). 

The amplitude of such a process, calculated in the fourth order in tunneling 
matrix elements, diverges logarithmically when the energy lo of an incoming electron 
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approaches 0. Since w ~ T, the logarithmic singularity in the transmission amplitude 
translates into a dramatic enhancement of the conductance G across the dot at low 
temperatures: G may reach values as high as the quantum limit 2e^//i |31[ \62\ . This 
conductance enhancement is not really a surprise. Indeed, in the spin-flip co-tunneling 
process a quantum dot with odd TV behaves as S = 1/2 magnetic impurity embedded 
into a tunneling barrier separating two massive conductors . It is known '341 since 
mid-60's that the presence of such impurities leads to zero-bias anomalies in tunneling 
conductance |SS| j which are adequately explained jSEl EZI in the context of the Kondo 
effect 0. 

At energies well below the threshold A ^ &E for intradot excitations the 
transitions within the (25' -I- l)-fold degenerate ground state manifold of a dot can 
be conveniently described by a spin operator S. The form of the effective Hamiltonian 
describing the interaction of the dot with conduction electrons in the leads is then 
dictated by SU{2) symmetry ||, 

^eft = 4clfesC„fe, + J^c^'iSa'c. ' S) (5.1) 

aks aa' 

with Saa' = J2kk'ss' ^aksi'^ ss' /2) c^,^,^,. The sum over k in Eq. H5.1|l is restricted to 
l^fel < A. The exchange amplitudes J^^, form a 2 x 2 Hermitian matrix J. The matrix 
has two real eigenvalues, the exchange constants Ji and J2 (hereafter we assume that 
Ji > J2) By an appropriate rotation in the R — L space the Hamiltonian (|5.2() can 
then be brought into the form 

fks 7 

Here the operators ip-y are certain linear combinations of the original operators cr^l 
describing electrons in the leads, and 

S7 = E '('Iks^^-yk's' 
kk' ss' 

is local spin density of itinerant electrons in the "channel" 7 = 1,2. 

The symmetry alone is not sufhcient to determine the exchange constants X,,; 
their evaluation must rely upon a microscopic model. Here we breifly outline the 
derivation |^ |2H1 of Eq. IjS.lfl for a generic model of a quantum dot system discussed 
in Sec. |21 above. For simplicity, we will assume that the gate voltage A'o is tuned to 
the middle of the Coulomb blockade valley. The tunneling H2.17|l mixes the state with 
N = Nq electrons on the dot with states having N ± 1 electrons. The electrostatic 
energies of these states are high (~ Eq), hence the transitions iV ^ TV ± 1 are 
virtual, and can be taken into account perturbatively in the second order in tunneling 
amplitudes [5^ . 

For Hamiltonian (|2.13() the occupations of single-particle energy levels are good 
quantum numbers. Therefore, the amplitude Jaa' can be written as a sum of partial 
amplitudes, 

Jaa' = Jga' ■ (5-3) 

n 

II In writing Eq. 15. li wo omitted the potential scattering terms associated with the usual elastic 
co-tunneling. This approximation is well justified when the conductances of the dot-lead junctions 
are small, Ga <IC e^/h, in which case Gei is also very small, see Eq. 14.51 
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Each term in the sum here corresponds to a process during which an electron or a hole 
is created virtually on level n in the dot, cf. Eq. (|4.4|) . For ^ /h and Es <C 5E 
the main contribution to the sum in H5.3|l comes from singly-occupied energy levels 
in the dot. A dot with spin S has 25 such levels near the Fermi level (hereafter we 
assign indexes n = — S*, . . . , n = S* to these levels), each carrying a spin S/25, and 
contributing 

Jaa' — -E^^*an^a'm Ki — '^/S, \n\ < S (5.4) 

Ec 

to the exchange amplitude in 1)5. 1() . This yields 

Jaa' ~ Jaa'- (5-5) 

|ri|<S 

It follows from Equations 1)5. 3|) and 15. 4|) that 



1 



trJ=— >^A„(|iL| + |4J)- (5.6) 



By restricting the sum over n here to |n| < S, as in H5.5|) . and taking into account 
that all A„ in H5.4|l are positive, we find Ji + J2 > 0. Similarly, from 

i = ^ E ^^>^n\VU, V^n = det ( \ (5.7) 

and Equations (|5.4I) and H5.5|l it follows that Ji J2 > for S* > 1/2. Indeed, in this case 
the sum in H5.7|l contains at least one contribution with m ^ n\ all such contributions 
are positive. Thus, both exchange constants Ji_2 > if the dot's spin S exceeds 
1/2 [23. The pecularities of the Kondo effect in quantum dots with large spin are 
discussed in Sec. [7|below. 

Here we concentrate on the most common situation of S* = 1/2 on the dot 
considered in detail in Sec. El The ground state of such dot has only one singly- 
occupied energy level [n = 0), so that det J w 0, see (|5.5|l and l|5.7(l . Accordingly, one 
of the exchange constants vanishes, 

J2 ~ 0, (5.8) 

while the remaining one, Ji = tr J, is positive. Equation (|5.8|) resulted, of course, 
from the approximation made in 1)5. 5|l . For the model 1)2. 13(1 the leading correction 
to ((5.5(1 originates in the co-tunneling processes with an intermediate state containing 
an extra electron (or an extra hole) on one of the empty (doubly-occupied) levels. 
Such contribution arises because the spin on the level n is not conserved by the 
Hamiltonian 1(2.13(1 . unlike the corresponding occupation number. Straightforward 
calculation 38. yields the partial amplitude in the form of Eq. ((5.4(1 , but with 

2EcEs 

{E^+W\f 

Unless the tunneling amplitudes tao to the only singly-occupied level in the dot 
are anomalously small, the corresponding correction 

5Jaa' = E J-^a' (5-9) 

to the exchange amplitude 1(5.5(1 is small. 



, I, n2 ' n^Q. 



OiCt' 



J OL 
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see Eq. (|2.14ll . To obtain this estimate, we assumed that ah tunnehng amphtudes 
tan are of the same order of magnitude, and replaced the sum over n in (|5.9|l by an 
integral. A similar estimate yields the leading contribution to det J, 

det J « ^ AoA„|I?o'J ^ (tr i)', 

C n 

or 

J2/J1 - -Es/5E. (5.10) 

According to (|5.10|l . the exchange constant J2 is negative 001 1 and its absolute 
value is small compared to Ji. Hence, 1)5.811 is indeed an excellent approximation for 
large chaotic dots with spin S = 1/2 as long as the intradot exchange interaction 
remains weak, Es ^ SE%. Note that corrections to the universal Hamiltonian ()2.13|l 
also result in finite values of both exchange constants, | J2I ~ JiN~^^'^, and become 
important for small dots with < 10 ,32^. Although this may significantly affect 
the conductance across the system in the weak coupling regime T > Tk, it does not 
lead to qualitative changes in the results for S — 1/2 on the dot, as the channel 
with smaller exchange constant decouples at low energies 021 j see also Seed below. 
With this caveat, we adopt approximation (|5.8|l in our description of the Kondo effect 
in quantum dots with spin S = 1/2. Accordingly, the effective Hamiltonian of the 
system (|5.2|l assumes the "block-diagonal" form 

H,tf = Hi+H2 (5.11) 

Hi=J2^kAs^ik. + Ji^i-S) (5.12) 

ks 
ks 

with J = tr J > 0. 



6. Kondo regime in transport through a quantum dot 

To get an idea about the physics of the Kondo model (see for recent reviews) , let 
us first replace the fermion field operator si in Eq. H5.12|l by a single-particle spin-1/2 
operator Si. The ground state of the resulting Hamiltonian of two spins 

if = J(Si -S) 

is obviously a singlet. The excited state (a triplet) is separated from the ground state 
by the energy gap Ji . This separation can be interpreted as the binding energy of the 
singlet. Unlike Si in this simple example, the operator si in H5.12|l is merely a spin 
density of the conduction electrons at the site of the "magnetic impurity" . Because 
conduction electrons are freely moving in space, it is hard for the impurity to "capture" 
an electron and form a singlet. Yet, even a weak local exchange interaction suffices to 
form a singlet ground state |44l I45| . However, the characteristic energy (an analogue 
of the binding energy) for this singlet is given not by the exchange constant J, but by 
the so-called Kondo temperature 

- Aexp(-l/i/J). (6.1) 

^ Equation 15. 8i holds identically for the Anderson impurity model 1371 frequently employed to study 
transport through quantum dots |81II41I . In that model a quantum dot is described by a single energy 
level, which formally corresponds to the infinite level spacing limit 5E — > 00 of the Hamiltonian I2.13i . 
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Using A ~ 6E and Equations H5.6|l and H2.19|l . one obtains from H6.1|l the estimate 
(5E\ 1 e^/h Ec 

\t^)^-j^ gTTg-.Je- (^-2) 

Since Gq ^ e'^/h and Eq 3> (J-E, the r.h.s. of (|6.2|l is a product of two large 
parameters. Therefore, the Kondo temperature Tk is small compared to the mean 
level spacing, 

Tk < SE. (6.3) 

It is this separation of the energy scales that justifies the use of the effective low- 
energy Hamiltonian H5.1|l . (|5.2|) for the description of the Kondo effect in a quantum 
dot system. Inequality H6.3|l remains valid even if the conductances of the dot-leads 
junctions Ga are of the order of 2e^//i. However, in this case the estimate (|6.2ll is no 
longer apphcable |i46j . 

In our model, see Equations H5.11(l - I|5.13(l . one of the channels (V'2) of conduction 
electrons completely decouples from the dot, while the ■i/'i-particles are described by 
the standard single-channel antiferromagnetic Kondo model j^lllSI- Therefore, the 
thermodynamic properties of a quantum dot in the Kondo regime are identical to those 
of the conventional Kondo problem for a single magnetic impurity in a bulk metal; 
thermodynamics of the latter model is fully studied |47j . However, all the experiments 
addressing the Kondo effect in quantum dots test their transport properties rather 
then thermodynamics. The electron current operator is not diagonal in the {'ipi,'ip2) 
representation, and the contributions of these two sub-systems to the conductance are 
not additive. Below we relate the linear conductance and, in some special case, the 
non-linear differential conductance as well, to the t-matrix of the conventional Kondo 
problem. 

6.1. Linear response 

The linear conductance can be calculated from the Kubo formula 

where the current operator / is given by 

N^^Y.'^l.^^o.Us (6-5) 

ks 

Here Na is the operator of the total number of electrons in the lead a. Evaluation of the 
linear conductance proceeds similarly to the calculation of the impurity contribution 
to the resistivity of dilute magnetic alloys (see, e.g., :48 ). In order to take the full 
advantage of the decomposition H5.11|l - H5.13|l . we rewrite / in terms of the operators 
4>i,2- These operators are related to the original operators cji^l representing the 
electrons in the right and left leads via 

V'lfcs \ ^ ( cos 6*0 sin 6*0 \ / CRks 
i^2ks J V -sin 6*0 cos 6*0 / V CLks 



(6.6) 



The rotation matrix here is the same one that diagonalizes matrix J of the exchange 
amplitudes in H5.1|) : the rotation angle satisfies the equation tan^o = |^io/^iio|- 
With the help of Eq. (|6.6|l we obtain 

Nr~Nl = cos(20o)(iVi - N2) - sin(20o) ^(^L.V'afcs + H-C.) (6.7) 



ks 
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The current operator / entering the Kubo formula (|6.4|) is to be calculated with 
the equilibrium Hamiltonian (|5.11|) - (|5.13|) . Since both Ni and N2 commute with 
H^fi, the first term in H().7|l makes no contribution to /. When the second term 
in Ht).7|l is substituted into (|6.5|) and then into the Kubo formula H6.4|) . the result, 
after integration by parts, can be expressed via 2-particle correlation functions such 
as (V'l (i)V'2(*)V'2(0)V'i(0)) (see Appendix B of 0^1 for further details about this 
calculation). Due to the block-diagonal structure of -ffgff, see (|5.11|l . these correlation 
functions factorize into products of the single-particle correlation functions describing 
the (free) ?/'2 -particles and the (interacting) T/ii-particles. The result of the evaluation 
of the Kubo formula can then be written as 

G = GoJ (^-^] i^[-7rz.ImT,(u;)]. (6.8) 

Here 



Go = ^ sin2(2^^o) = (6 9) 

/(cj) is the Fermi function, and Ts{u!) is the t-matrix for the Kondo model (|5.12|) . 
The t-matrix is related to the exact retarded Green function of the '0i-particles in the 
conventional way, 

GksM^) = Gli^) + GI{lo)Ts{uj)GI,{uj), Gl = {u - & + *0)-i. 

Here Gks,k's{^) is the Fourier transform of Gks,k's{t) ~ ^*^(i)({''/'ifcs(^)i ^ik's})' where 
(. . .) stands for the thermodynamic averaging with Hamiltonian H5.12|l . In writing 
Eq. H6.8|l we took into account the conservation of the total spin (which implies that 
Gks,k's' = Sss'Gks,k'sy and that the interaction in H5.12|l is local (which in turn means 
that the t-matrix is independent of k and k'). 

6. 2. Weak coupling regime: Tk <C T <C 5E 

When the exchange term in the Hamiltonian H5.12|l is treated perturbatively, the main 
contribution to the t-matrix comes from the transitions of the type j5()| 

\ks,(j) ~^\k's\a') . (6.10) 

Here state |fcs, a) has an extra electron with spin s in the orbital state k whereas the 
dot is in the spin state a. By SU (2) symmetry, the amplitude of the transition (|6.1U|I 
satisfies 

A\k's',(j')^\ks,(T) = Ak'k ^ (Cs's • CTa'a) (6-11) 

The transition (|6.10(l is elastic in the sense that the number of quasiparticles in the 
final state of the transition is the same as that in the initial state (in other words, 
the transition l|6.10|) is not accompanied by the production of electron-hole pairs). 
Therefore, the imaginary part of the t-matrix can be calculated with the help of the 
optical theorem [HJ, which yields 



jEEk^4'cv,.')-ife«,-> ■ (6.12) 



— ■KV Im T, = 

2 

(7 S'a' 

The factor 1/2 here accounts for the probability to have spin a on the dot in the initial 
state of the transition. Substitution of the tunneling amplitude in the form (|6.11|l into 
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Eq. Ht).12|l . and summation over the spin indexes with the help of the identity H2.7|l 
result in 

-T,vlTnTs=^-^v^\Al,^\. (6.13) 

The amplitude Ak'k in Equations Hti.ll|l and ()6.13f) depends only on the difference 
of energies lo = £,k' - Cfe: 

Ak'k = A{uj). 

In the leading (first) order in J one readily obtains A^^^ = J, independently of u. 
However, as discovered by Kondo the second-order contribution A^'^'> not only 
depends on but is logarithmically divergent as cj ^ 0: 

A(2)(tj) \n\l\luo\. 

Here A is the high-energy cutoff in the Hamiltonian H5.12|l . It turns out 'SO' that 
similar logarithmically divergent contributions appear in all orders of perturbation 
theory, 

i/A(")(lj) = (!.J)"[ln|A/c^|]""\ 
resulting in a geometric series 

CO OO J 

Ti=l ri=0 ' ' ' 

With the help of the definition of the Kondo temperature l|6.1l) . this can be written as 

vA{^) = / . (6.14) 
In \ujITk\ 

Substitution of H6.14|l into Eq. H6.13|l and then into Eq. 16.811 , and evaluation of the 
integral over with logarithmic accuracy yield for the conductance across the dot 

3^716 

•\j,^{TITk) 

Equation H6.15|l is the leading term of the asymptotic expansion in powers 
of l/ln(r/r/f), and represents the conductance in the leading logarithmic 
approximation. 

Equation (|6.15() resulted from summing up the most-diverging contributions in 
all orders of perturbation theory. It is instructive to re-derive it now in the framework 
of renormalization group (52) . The idea of this approach rests on the observation that 
the electronic states that give a significant contribution to observable quantities, such 
as conductance, are states within an interval of energies of the width w ~ T about 
the Fermi level, see Eq. (|6.8|l . At temperatures of the order of Tk, when the Kondo 
effect becomes important, this interval is narrow compared to the width of the band 
D = A in which the Hamiltonian H5.12|l is defined. 

Consider a narrow strip of energies of the width 5D ^ D near the edge of the 
band. Any transition H6.10|l between a state near the Fermi level and one of the states 
in the strip is associated with high (~ A) energy deficit, and, therefore, can only occur 
virtually. Obviously, virtual transitions via each of the states in the strip result in the 
second-order correction ^ / D to the amplitude A{uj) of the transition between 
states in the vicinity of the Fermi level. Since the strip contains vSD electronic states, 
the total correction is 52 

5 A ~ vj'^5D/D. 



G = Go ^ ^J,^ , , T»Tif. (6.15) 
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This correction can be accounted for by modifying the exchange constant in the 
effective Hamiltonian -ffeff which is defined for states within a narrower energy band 
of the width D - 5D 

HcH^Y.^ki^is'^ks + JD-SDis^-S), \^k\<D-SD, (6.16) 

ks 

Jd-sd = Jd + i^Je)-^- (6-17) 

Here Ju is the exchange constant in the original Hamihonian. Note that the -ffcff 
has the same form as Eq. (|5.12|) . This is not merely a conjecture, but can be shown 
rigorously 021 EHl- 

The reduction of the bandwidth can be considered to be a result of a unitary 
transformation that decouples the states near the band edges from the rest of the 
band. In principle, any such transformation should also affect the operators that 
describe the observable quantities. Fortunately, this is not the case for the problem 
at hand. Indeed, the angle 6*0 in Eq. l|6.6|l is not modified by the transformation. 
Therefore, the current operator and the expression for the conductance H6.8I) retain 
their form. 

Successive reductions of D by small steps SD can be viewed as a continuous 
process during which the initial Hamiltonian (|5.12(l with _D = A is transformed to an 
effective Hamiltonian of the same form that acts within the band of the reduced width 
D A. It follows from H6.17|l that the dependence of the effective exchange constant 
on D is described by the differential equation El] 

'^ = ^Jl, C = ln(A/i?). (6.18) 

With the help of Eq. H6.1|l . the solution of the RG equation (|6.18|l subject to the initial 
condition — J can be cast into the form 



Hd/Tk) ■ 

The renormalization described by Eq. H6.18|l can be continued until the bandwidth D 
becomes of the order of the typical energy |cj| ~ T for real transitions. After this limit 
has been reached, the transition amplitude A{lu) is calculated in lowest (first) order 
of perturbation theory in the effective exchange constant (higher order contributions 
are negligibly small for D ^ lu ), 

ln|w/TA'| 

Using now Equations 16.13|l and (|6.8|l . we recover Eq. I|6.15|l . 



6.3. Strong coupling regime: T <C Tji 

As temperature approaches Tk, the leading logarithmic approximation result (|6.15|l 
diverges. This divergence signals the failure of the approximation. Indeed, we are 
considering a model with single-mode junctions between the dot and the leads. The 
maximal possible conductance in this case is 2e^//i. To obtain a more precise bound, 
we discuss in this section the conductance in the strong coupling regime T <^ Tk- 

We start with the zero-temperature limit T — 0. As discussed above, the ground 
state of the Kondo model H5.12|l is a singlet 01], and, accordingly, is not degenerate. 
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Therefore, the t-matrix of the conduction electrons interacting with the locahzed spin 
is completely characterized by the scattering phase shifts 5^ for electrons with spin 
s at the Fermi level. The t-matrix is then given by the standard scattering theory 
expression |SJ 

-^yTs[Q)^^{^s-l). (6.19) 

where §s is the scattering matrix for electrons with spin s, which for a single channel 
case reduces to its eignevalue. Substitution of I|fi.l9|l into Eq. Hfi.8|l yields 

G(0) = Goi^sin2,5, (6.20) 

s 

for the conductance, see Eq. H6.8|l . The phase shifts in (|6.19(l . (|6.2()|l are obviously 
defined only modvr (that is, 5s and <5s + tt are equivalent). This ambiguity can be 
removed if we set to zero the values of the phase shifts at J = in Eq. (|5.12(l . 

In order to find the two phase shifts 5s, we need two independent relations. 
The first one follows from the invariance of the Kondo Hamiltonian (|5.12() under the 
particlc-holc transformation ■0^3 ^"^-k -s (here s = ±1 for spin-up/down electrons). 
The particle-hole symmetry implies the relation for the t-matrix 

Tsiiu) = ~Tl,{-uj), (6.21) 

valid at all w and T. In view of Eq. H6.19|l . it translates into the relation for the phase 
shifts at the Fermi level {u = 0) j54| . 

+ -^i = 0. (6.22) 

The second relation follows from the requirement that the ground state of the 
Hamiltonian (|5.12(l is a singlet In the absence of exchange ( J = 0) and at T = 0, 
an infinitcsimally weak {B — > -1-0) magnetic field acting on the dot's spin, 

Hb = BS^, (6.23) 

would polarize it; here B is the Zeeman energy. Since free electron gas has zero spin in 
the ground state, the total spin in a very large but finite region of space V surrounding 
the dot coincides with the spin of the dot, {S^)j=q = —1/2. If the exchange with 
electron gas is now turned on, J > 0, a very weak field will not prevent the formation 
of a singlet ground state. In this state, the total spin within V is zero. Such change of 
the spin is possible if the numbers of spin-up and spin-down electrons within V have 
changed to compensate for the dot's spin: (5iV| — 5Ni — 1. By the Friedel sum rule, 
5Ns are related to the scattering phase shifts at the Fermi level, 5Ns — 5s /tt, which 
gives 

5^-Si^ TT. (6.24) 

Combining and we find \5s\ = 7r/2. Equation if^?^ then yields for 

zero-temperature and zero-field conductance across the dot 

G(0) = Go. (6.25) 

Thus, the grow of the conductance with lowering the temperature is limited only by 
the value of Go- This value, see Eq. H6.9|l . depends only on the ratio of the tunneling 
amplitudes \tLo/tRo\. If \tLo\ = \tRo\, then the conductance at T = will reach the 
maximal value G — 2e'^/h allowed by quantum mechanics 13 Ij . 

The maximal conductance, Eq. H6.25|l . is reached when a a singlet state is 
formed by the itinerant electrons interacting with the local spin, as described by 



Kondo effect in quantum dots 



18 



G 

Go 

y 

(6.27) 

Goc 
Gel 





1 (6.15) 


(4.5) 


5' = 

(4.2) 


1 
2 

(3.6) 


(3.1) 























Tk SE Tel 



Ec 



T 



Figure 4. Sketch of the temperature dependence of the conductance in the 
middle of the Coulomb blockade valley with S = 1/2 on the dot. The numbers in 
brackets refer to the corresponding equations in the text. 

the Kondo Hamiltonian (|5.12|1 . Perturbation of this singlet [S^ by a magnetic field 
B or temperature T leads to a decrease of the conductance. This decrease is small as 
long as B and T are small compared to the singlet "binding energy" Tk- The reader 
is referred to the original papers [Si] for the details. Here we only quote the result 0H1 
for the imaginary part of the t-matrix at and T small compared to the Kondo 
temperature Tk, 

~ T:iyImTs{u;) = 1 - ^"^^ ^jj^^ ■ (6.26) 



2T^ 



Substitution of into (|^ yields 



G — Go 



1 - {ttT/TkY 



T-^Tk. 



(6.27) 



Accordingly, corrections to the conductance are quadratic in temperature - a typical 
Fermi liquid result |54) . The weak-coupling (T 3> Tk) and the strong-coupling 
(T <ti Tk) asymptotes of the conductance have a strikingly different structure. 
Nevertheless, since the Kondo effect is a crossover phenomenon rather than a phase 
transition (43i .411 1451 14 7| . the dependence G{T) is a smooth and featureless 
function throughout the crossover region T ^ Tk- 

Finally, note that both Eq. I|6.15|l and Eq. H6.27|) have been obtained here for the 
particle-hole symmetric model (|5.12(l . This approximation is equivalent to neglecting 
the elastic co-tunneling contribution to the conductance Gei- The asymptotes (|6.15|l 
and (|6.27|) remain valid as long as Gei/Go <^ 1. The overall temperature 
dependence of the linear conductance in the middle of the Coulomb blockade valley is 
sketched in Fig. ^ 

6-4- Beyond linear response 

In order to study transport through a quantum dot away from equilibrium we add to 
the effective Hamiltonian (|5.11|) - (I5.13|) a term 

Hv = ^ {Nl - Nnj (6.28) 
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describing a finite voltage bias V applied between the left and right electrodes. Here 
we will evaluate the current across the dot at arbitrary V but under the simplifying 
assumption that the dot-lead junctions are strongly asymmetric: 

Gl ^ Gr. 

Under this condition the angle 6*0 in l|6.6(l is small, 6*0 ~ \tLo/tRa\ <C 1. Expanding 
Eq. (|6.7|l to linear order in we obtain 

Hv{0^) = ^(^2 - ^i) + ey^o E('/'lfc^^2fe. + H.c.) (6.29) 

ks 



The first term in the r.h.s. here can be interpreted as the voltage bias between 
the reservoirs of 1- and 2-particles, cf. Eq. H6.28|l . while the second term has an 
appearance of fc-conserving tunneling with very small (proportional to ^ 1) 
tunneling amplitude. 

Similar to Eq. H6.29|l . the current operator /, see (|6.5|) . splits naturally into two 
parts, 

/ = /o + 51, 

ks 

ks 

It turns out that 61 does not contribute to the average current in the leading (second) 
order in 9o The remaining contribution / = (/q) corresponds to tunneling current 
between two bulk reservoirs containing 1- and 2-particles. Its evaluation yields 

^^Golj2^-n,.ImTs{eV)] (6.30) 

s 

for the differential conductance across the dot at zero temperature. Here Gq coincides 
with the small 6'o-hmit of Eq. (|^^ . Using now Equations and 
we obtain 

(6.31) 

eV:^TK 



1 dl \ 2 \TkJ 
G^dV ^ I 3^Vl6 



ln\eV/TK) ' 

Thus, a large voltage bias has qualitatively the same destructive effect on the 
Kondo physics as the temperature does. The result (|6.31|) remains valid as long as 
T <C eV <C SE. If temperature exceeds the bias, T ^ eV, the differential conductance 
coincides with the linear conductance, see Equations Ijfi.lSII and (|6.27|) above. 



7. Kondo effect in quantum dots with large spin 

If the dot's spin exceeds 1/2 |56l I57L then, as discussed in Sec. HI above, both 
exchange constants in the effective Hamiltonian (|5.2|l are finite and positive. This 
turns out to have a dramatic effect on the dependence of the conductance in the Kondo 
regime on temperature T and on Zeeman energy B. Unlike the case of 5 = 1/2 on the 
dot, see Fig.^ now the dependences on T and B are nonmonotonic: initial increase of 
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G follows by a drop when the temperature is lowered |^ EH] at i? = 0; the variation 
of G with B at T = is similarly nonmonotonic. 

The origin of this peculiar behavior is easier to undersand by considering the in- 
dependence of the zero-temperature conductance |211- We assume that the magnetic 
field H\\ is applied in the plane of the dot. Such field leads to the Zeeman splitting B 
of the spin states of the dot, see Eq. H6.23(l . but barely affects the orbital motion of 
electrons. 

At any finite B the ground state of the system is not degenerate. Therefore, the 
linear conductance at T = can be calculated from the Landauer formula 

G=^^Y.Krl\. (7.1) 

which relates G to the amplitude of scattering §s;_rl of an electron with spin s from 
lead L to lead R. The amplitudes §s;Qa' form a 2 x 2 scattering matrix S^. In the 
basis of "channels", see Eq. 1)5.2(1 . this matrix is obviously diagonal, and its eigenvalues 
exp(2i(5-ys) are related to the scattering phase shifts The scattering matrix in the 
original {R — L) basis is obtained from 

^ C/tdiag {e^**-'^ } U, 

where [/ is a matrix of a rotation by an angle Oq, see Eq. (|5.t)|) . The Landauer 
formula (|7.1() then yields 

G = Goi5]sin2(5i, -<52s), Go = ^sin2(20o), (7.2) 

s 

which generalizes the single-channel expression H6.2U|I . 

Equation (|7.2() can be further simplified for a particle-hole symmetric model (|5.2|) . 
Indeed, in this case the phase shifts satisfy S^i+S^i — 0, cf. Eq. 1(6.22)1 . which suggests 
a representation 

Substitution into 1(7.2(1 then results in 

G = Go sin^ {6i - 62) . (7.3) 




Figure 5. Dependence of the scattering phase shifts at the Fermi level on the 
magnetic field for S > 1/2 (left panel) and S = 1/2 (right panel). 

If spin on the dot S exceeds 1/2, then both channels of itinerant electrons 
participate in the screening of the dot's spin Accordingly, in the limit B ^ both 
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phase shifts 5^ approach the unitary hmit value 7r/2, see Fig.|Sl However, the increase 
of the phase shifts on lowering the field is characterized by two different energy scales. 
These scales, the Kondo temperatures Ti and T2, are related to the corresponding 
exchange constants in the effective Hamiltonian 1)5.2(1 . 




so that Ti > T2 for Ji > J2- It is then obvious from Eq. (|7.3() that the conductance 
across the dot is small both at weak {B <C T2) and strong {B 3> Tf) fields, but 
may become large (~ Go) at intermediate fields T2 ^ S ^ Ti, see Fig. El In 
other words, the dependence of zero-temperature conductance on the magnetic field 
is nonmonotonic. 




Figure 6. Sketch of the magnetic field dependence of the Kondo contribution 
to the linear conductance at zero temperature. The conductance as function of 
temperature exhibits a similar nonmonotonic dependence. 



This nonmonotonic dependence is in sharp contrast with the monotonic increase 
of the conductance with lowering the field when S = 1/2. Indeed, in the latter case it 
is the channel whose coupling to the dot is the strongest that screens the dot's spin, 
while the remaining channel decouples at low energies see Fig. |3 

The dependence of the conductance on temperature G{T) is very similar to 
G{B)^ . For example, for 5 = 1 one obtains 

2 



G/Gq — 



1 



1 



Y 



Ti 
1 



T2 



1 



T<T2 



(7.4) 



ln(T/Ti) ln(T/T2)_ 

VTfT^. The value of G 

max 

can be found analytically for Ti ^ T2. For S = 1 the result reads [231 

3^2 - 



The conductance reaches its maximal value Gmax a-t T 



G„ 



Go 



1 



ln^(Ti/T2) 



(7.5) 



Consider now a Coulomb blockade valley with N — even electrons and spin S = \ 
on the dot. In a typical situation, the dot's spin in two neighboring valleys (with iV± 1 
electrons) is 1/2. Under the conditions of applicability of the approximation (|5.5|) . 
there is a single non-zero exchange constant Jn±i for each of these valleys. If 
the Kondo temperatures Tx are the same for both valleys with S = 1/2, then 
Jn+1 = Jn-1 = Jodd- Each of the two singly-occupied energy levels in the valley 

+ Note , how ever, that (Vij (i)V'2 W'/'l (O)V'i (0)> 7^ {i;l{t)ip^{0)){^lj^{t)ipl{0)) at finite T. Therefore, 
unlike I6.2UI . Eq. 16.81 does not allow for a simple generalization to the two-channel case. 



Kondo effect in quantum dots 



22 



with S" = 1 is also singly-occupied in one of the two neighboring valleys. It then 
follows from Equations H5.4|) - (|5.t)|) that the exchange constants Ji 2 for S* = 1 satisfy 

1 

+ = 2 [Jn+I + Jn-i) = Jodd- 

Since both Ji and J2 are positive, this immediately implies that Ji^2 < Jodd- 
Accordingly, both Kondo temperatures Ti^2 are expected to be smaller than Tk in 
the nearby valleys with S = 1/2. 

This consideration, however, is not applicable when the dot is tuned to the vicinity 
of the singlet-triplet transition in its ground state [J!, 'W, "37", , i.e. when the energy 
gap A between the triplet ground state and the singlet excited state of an isolated dot 
is small compared to the mean level spacing SE. In this case the exchange constants 
in the effective Hamiltonian H5.2|) should account for additional renormalization that 
the system's parameters undergo when the high-energy cutoff (the bandwidth of the 
effective Hamiltonian) D is reduced from D ^ SE down to 13 ^ A <C 6E j^H], see 
also The renormalization enhances the exchange constants Ji.2. If the ratio 

A/SE is sufficiently small, then the Kondo temperatures Ti^2 for S ^ 1 may become 
of the same order [^^1 IHHI , or even significantly exceed [3 IHl EI] the corresponding 
scale Tk for S* = 1/2. 

In GaAs-based lateral quantum dot systems the value of A can be controlled by 
a magnetic field H± applied perpendicular to the plane of the dot |57) . Because of 
the smallness of the effective mass m*, even a weak field H± has a strong orbital 
effect. At the same time, smallness of the quasiparticle g-factor in GaAs ensures that 
the corresponding Zeeman splitting remains small [HI. Theory of the Kondo effect 
in lateral quantum dots in the vicinity of the singlet-triplet transition was developed 
in [HT], see also |^ . 



8. Discussion 



In the simplest form of the Kondo effect considered in this review, a quantum 
dot behaves essentially as an artificial "magnetic impurity" with spin S, coupled 
via exchange interaction to two conducting leads. The details of the temperature 
dependence G{T) of the linear conductance across a lateral quantum dot depend on 
the dot's spin S. In the most common case of S" = 1/2 the conductance in the Kondo 
regime monotonically increases with the decrease of temperature, potentially up to 
the quantum limit 2e^//i. Qualitatively (although not quantitatively), this increase 
can be understood from the Anderson impurity model in which the dot is described 
by a single energy level. On the contrary, when spin on the dot exceeds 1/2, the 
evolution of the conductance proceeds in two stages: the conductance first raises, and 
then drops again when the temperature is lowered. 

In a typical experiment [3], one measures the dependence of the differential 
conductance on temperature T, Zeeman energy B, and dc voltage bias V. When 
one of these parameters is much larger than the other two, and is also large compared 
to the Kondo temperature Tk, the differential conductance exhibits a logarithmic 
dependence 

ma.x{T,B,eV}^^^ 



1 dl 



In ^^^J ' (8-1) 

characteristic for the weak coupling regime of the Kondo system. Consider now a 
zero-temperature transport through a quantum dot with S = 1/2 in the presence of a 
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strong field B ^ T^. In accordance with H8.1|) . the differential conductance is small 
compared to Gq both for eV ^ B and for eV ^ B. However, the calculation in the 
third order of perturbation theory in the exchange constant yields a contribution that 
diverges logarithmically at eV = B The divergence appears because at eV — B 
the scattered electron has just the right amount of energy to allow for a real transition 
with a flip of spin. However, the full development of resonance is inhibited by a finite 
lifetime of the excited spin state of the dot |41[ Ift3| . As a result, the peak in the 
differential conductance at eV ^ B \s broader and lower 01] then the corresponding 
peak at zero bias in the absence of the field. Even though for eV ~ i? ^ Tr- the 
system is clearly in the weak coupling regime, a resummation of the perturbation 
series turns out to be a very difficult task, and the detailed shape of the peak is still 
unknown. This problem remains to be a subject of active research, see e.g. |()4| and 
references therein. 

One encounters similar difhculties in studies of the effect of a weak ac excitation 
of frequency > Tk applied to the gate electrode [SS] on transport across the dot. 
In close analogy with the usual photon-assisted tunneling such perturbation is 
expected to result in the formation of satellites [HZ| at eV = nMl (here n is an 
integer) to the zero-bias peak in the differential conductance. Again, the formation 
of the satellite peaks and the survival of the zero-bias peak in the presence of the ac 
excitation are limited by the finite lifetime effects |68| . 

The spin degeneracy is not the only possible source of the Kondo effect in quantum 
dots. Consider, for example, a large dot connected by a single-mode junction to a 
conducting lead and tuned to the vicinity of the Coulomb blockade peak |22|. If one 
neglects the finite level spacing in the dot, then the two almost degenerate charge 
state of the dot can be labeled by a pseudospin, while real spin plays the part of the 
channel index |22[ BU) . This setup turns out to be a robust realization |22[ IHU) of 
the symmetric (i.e. having equal exchange constants) two-channel S* = 1/2 Kondo 
model 021 ■ The model results in a peculiar temperature dependence of the observable 
quantities, which at low temperatures follow power laws with manifestly non-Fermi- 
liquid fractional powers. 

It should be emphasized that in the usual geometry consisting of two leads 
attached to a small* Coulomb-blockaded quantum dot with S = 1/2, only the 
conventional Fermi-liquid behavior can be observed at low temperatures. Indeed, 
in this case the two exchange constants in the effective exchange Hamiltonian H5.2|l 
are vastly different, and their ratio is not tunable by conventional means, see the 
discussion in Sec. [3 above. A way around this difficulty was proposed recently in |7Uj . 
The key idea is to replace one of the leads in the standard configuration by a very 
large quantum dot, characterized by a level spacing SE' and a charging energy E'q. 
At T ^ 6E' , particle- hole excitations within this dot are allowed, and electrons in it 
participate in the screening of the smaller dot's spin. At the same time, as long as 
T <C E'(j, the number of electrons in the large dot is fixed. Therefore, the large dot 
provides for a separate screening channel which does not mix with that supplied by 
the remaining lead. In this system, the two exchange constants are controlled by the 
conductances of the dot-lead and dot-dot junctions. A strategy for tuning the device 
parameters to the critical point characterized by the two-channel Kondo physics is 
discussed in |7T] . 

Finally, we should mention that the description based on the universal 

* i.e. with appreciable level spacing 
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Hamiltonian H2.13|l is not applicable to large quantum dots subjected to a quantizing 
magnetic field H±^ [721 HHI' Such field changes drastically the way the screening 
occurs in a confined droplet of a two-dimensional electron gas |74|. The droplet is 
divided into alternating domains containing compressible and incompressible electron 
liquids. In the metal-like compressible regions, the screening is almost perfect. On 
the contrary, the incompressible regions behave very much like insulators. In the 
case of lateral quantum dots, a large compressible domain may be formed near the 
center of the dot. This domain is surrounded by a narrow incompressible region 
separating it from another compressible ring-shaped domain formed along the edges 
of the dot |75| . This system can be viewed as two concentric capacitively coupled 
quantum "dots" - the core dot and the edge dot [73 CS]. When the leads are attached 
to the edge dot, the measured conductance is sensitive to its spin state: if the number 
of electrons in the edge dot is odd, then the conductance becomes large due to the 
Kondo effect [72 ■ Changing the field causes redistribution of electrons between the 
core and the edge, resulting in a striking checkerboard-like pattern of high- and low- 
conductance regions |72[ I73|. This behavior persists as long as the Zeeman energy 
remains small compared to the Kondo temperature. Note that compressible regions 
are also formed around an antidot - a potential hill in a two-dimensional electron gas 
in the quantum Hall regime [JHl- Both Coulomb blockade oscillations and Kondo- like 
behavior were observed in these systems too (77) . 

9. Summary 

Kondo effect arises whenever a coupling to a Fermi gas induces transitions within 
otherwise degenerate ground state multiplet of an interacting system. Both coupling to 
a Fermi gas and interactions are naturally present in a nanoscale transport experiment. 
At the same time, many nanostructures can be easily tuned to the vicinity of a 
degeneracy point. This is why the Kondo effect in its various forms often influences 
the low temperature transport in meso- and nanoscale systems. 

In this article we reviewed the theory of the Kondo effect in transport through 
quantum dots. A Coulomb-blockaded quantum dot behaves in many aspects as an 
artificial "magnetic impurity" coupled via exchange interaction to two conducting 
leads. Kondo effect in transport through such "impurity" manifests itself in the lifting 
of the Coulomb blockade at low temperatures, and, therefore, can be unambiguously 
identified. Quantum dot systems not only offer a direct access to transport properties 
of an artificial impurity, but also provide one with a broad arsenal of tools to tweak the 
impurity properties, unmatched in conventional systems. The characteristic energy 
scale for the intra-dot excitations is much smaller than the corresponding scale for 
natural magnetic impurities. This allows one to induce degeneracies in the ground 
state of a dot which are more exotic than just the spin degeneracy. This is only one 
out of many possible extensions of the simple model discussed in this review. 
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